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Abstract — Generalized comb filters (GCFs) are effi- 
cient anti-aliasing decimation filters with improved selec- 
tivity and quantization noise (QN) rejection performance 
around the so called folding bands with respect to 
classical comb filters. 

In this paper, we address the design of GCF filters 
by proposing an efficient partial polyphase architecture 
with the aim to reduce the data rate as much as possible 
after the EA A/D conversion. We propose a mathemat- 
ical framework in order to completely characterize the 
dependence of the frequency response of GCFs on the 
quantization of the multipliers embedded in the proposed 
filter architecture. This analysis paves the way to the 
design of multiplier-less decimation architectures. 

We also derive the impulse response of a sample 3rd 
order GCF filter used as a reference scheme throughout 
the paper. 

Index Terms — A/D converter, CIC-filters, comb, dec- 
imation, decimation filter, GCF, generalized comb filter, 
partial polyphase, polyphase, power-of-2, EA, sine filters. 



I. Introduction and problem formulation 

The design of computationally efficient decimation 
filters for E A A/D converters is a well-known research 
topic [1], [2], [3]. Given a base-band analog input 
signal x(t) with bandwidth [— f x , +f x ], a EA A/D 
converter of order B produces a digital signal x(nT s ) 
by sampling x(t) at rate f s = 2pf x ^> 2f x , whereby 
p is the so called oversampling ratio. The normalized 
maximum frequency contained in the input signal is 
defined as f c = fy- = and the digital signal x(nT s ) 
at the input of the first decimation filter has frequency 
components belonging to the range fd G [— f c , f c ] 
(where fd denotes the digital frequency). This setup is 
pictorially depicted in the reference architecture shown 
in Fig. Q] 

Owing to the condition p ^> 1, the decimation 
of an oversampled signal x(nT s ) is efficiently [4] 
accomplished by cascading two (or more) decimation 
stages, as highlighted in Fig. [U followed by a FIR filter 
which provides the required selectivity on the sampled 
signal x(tiTn) at baseband. The first decimation filter 
is usually an 7V-th order comb filter decimating by 
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Fig. 1. General architecture of a two-stage decimation chain for 
£ A A/D converters, along with a pictorial representation of the key 
frequency intervals to be carefully considered for the design of the 
first decimation stage. 



D [2], [3], [5], whereby the order N has to be greater 
or equal to B + 1 [1], [3]. 

The design of a multistage decimation filter for EA 
converters poses stringent constraints on the shape 
of the frequency response of the first decimation 
stage. Considering the scheme in Fig. Q] the fre- 
quency response H\{e^ u ) of the first decimation filter 
must attenuate the QN falling inside the so called 
folding bands, i.e. the frequency ranges defined as 
[s - /c! 5 + fc] with k = 1, Lf J if D is even, 
and k — 1,..., L"^5^"J f° r D odd (for conciseness, 
the set of values assumed by k will be denoted as 
Kk throughout the paper). The reason is that the 
EA QN falling inside these frequency bands, will 
fold down to baseband because of the sampling rate 
reduction by D in the first decimation stage, irremedia- 
bly affecting the signal resolution after the multistage 
decimation chain. This issue is especially important 
for the first decimation stage, since the QN folding 
down to baseband has not been previously attenuated. 
Fig.[T]also shows the frequency range [0, f c ] where the 
useful signal bandwidth falls, along with the so called 
don't care bands, i.e. the frequency ranges whose QN 
will be rejected by the filters placed beyond the first 



decimation stage. 

In connection with the first decimation stage in the 
multistage architecture shown in Fig. [T| the required 
aliasing protection around the folding bands is usually 
guaranteed by a comb filter, which provides an in- 
herent antialising function by placing its zeros in the 
middle of each folding band. We recall that the transfer 
function of a A-th-order comb filter is defined as [2] : 
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where D is the desired decimation factor. 

With this background, let us provide a quick survey 
of the recent literature related to the problem addressed 
here. This survey is by no means exhaustive and is 
meant to simply provide a sampling of the literature 
in this fertile area. 

A 3rd order modified decimation sine filter was 
proposed in [6], and still further analyzed in [7], 
[8]. The class of comb filters was then generalized 
in [9], whereby the authors proposed an optimization 
framework for deriving the optimal zero rotations of 
GCFs for any filter order and decimation factor D. 

Other works somewhat related to the topic addressed 
in this paper are [10]-[15]. In [10] and [11] authors 
proposed computational efficient decimation filter ar- 
chitectures for implementing classical comb filters. 
In [12] authors proposed the use of decimation sharp- 
ened filters embedding comb filters, whereas in [13] 
authors addressed the design of a novel two-stage 
sharpened comb decimator. In [14], authors proposed 
novel decimation schemes for EA A/D converters 
based on Kaiser and Hamming sharpened filters, then 
generalized in [15] for higher order decimation filters. 

The main aim of this paper is to propose a flexible, 
yet effective, partial polyphase architecture for imple- 
menting the GCF filters proposed in the companion 
paper [9]. To this end, we first recall the z-transfer 
function of GCF filters for completeness, and, then, 
provide a mathematical formulation for deriving the 
impulse response of this class of decimation filters. 
The latter is needed for deriving the polyphase compo- 
nents of the proposed filters. In the second part of the 
paper, the focus is on the sensitivity of the frequency 
response of GCF filters due to the quantization of 
the multipliers embedded in the proposed architecture. 
We also analyze zero displacements in the z-transfer 
function of GCF filters for deducing useful hints at the 
basis of any practical implementation of such filters. 

The rest of the paper is organized as follows. In 
Section II, we briefly recall the transfer functions of 
GCF filters and highlight their main peculiarities with 
respect to classical comb filters. Section III presents an 
effective architecture, namely a partial polyphase de- 
composition, for implementing GCF decimation filters; 
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Fig. 2. Zero locations of the considered decimation filters within 
the fc-th folding band. Zeros are displayed in the z-plane and in the 
frequency axis in order to highlight their effect in both domains. 



the impulse response of a sample 3rd order GCF filter 
is also presented, and mathematically derived in the 
Appendix. In Section IV, we present a mathematical 
framework for evaluating the sensitivity of the fre- 
quency response of GCF filters to the approximations 
of the embedded multipliers. Zero displacements due 
to multiplier approximation is discussed in Section V, 
where we also draw general guidelines for the design 
of such filters. Finally, Section VI draws the conclu- 



II. The z-Transfer function of GCF filters 

In this section, we briefly recall the z-transfer func- 
tion of GCF filters proposed in [9]. The z-transfer 
function Hqcf n (z) of a A-th order GCF filter deci- 
mating by D is defined as: 

1 If I 

Hi{z) ■ rin=i Hi{z, otn+N), D even, N even 
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Hi(z) -J]^ H 2 {z,a n+N ), D even, N odd 



D odd (2) 



whereby the involved basic functions are defined as 
follows: 

hi(z) = nfj? n£Li (i - 2cos - °») * _1 + z ~ 2 ) > 

= nf=i ( X " ^" 1 e +J * 1 " iQ ") (l - z^e" 3 '*^'"' 

(3) 

with D M = ^ - 1, for D even, and D M = ^f^, for 
D odd, and 

H 2 (z, a n+N ) = 1-2 cos (it - a n+N ) z~ x + z~ 2 (4) 
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Terms H Q 



and H 0t0 d are appropriate 



normalization constants chosen in such a way as to 
have HccF N (z)\ z=1 = 1, and |_-J is the floor of the 
underlined number. 

Let us summarize the main peculiarities of GCF 
filters by comparing them to classical comb filters. 
GCFs are, as comb filters, linear-phase filters since 
they are constituted by two linear-phase basic filters, 
namely H\(z) and H-2.{z). 

An iV-th order comb filter (Q} decimating by D, 
places iV-th order zeros in the complex locations Zi — 
e- 7 "^ 1 , Vi = 1, . . . , D— 1, or, equivalently, in the digital 
frequencies f Zk = Mj, k G K^. On the other hand, 
an iV-th order GCF filter decimating by D places N 
pairs of conjugate complex zeros in the i-th folding 
band! with i = 1, . . . , Dm, as exemplified by the 
function H\(z) in (0. A pictorial representation of the 
zero locations of a GCF filter is given in Fig. [2] The 
behaviour of the function H2(z) is analogous to that 
of Hi(z) with the exception that its zeros are placed 
around the location z = — 1 in the z-complex plane, 
and, it holds only for D even. A convenient choice for 
a p is a p = q p 2irf c , with q p G [— 1,+1]: this solution 
is such that each pair of conjugate complex zeros 
falls inside the relative folding band guaranteeing the 
required selectivity in these frequency bands. By virtue 
of the zero distribution within the folding bands, GCF 
filters provide improved SA QN rejection capabilities 
with respect to classical comb filters of the same 
order. For completeness, Table H] shows the optimal 
zero rotations q p s found in [9] by minimizing the 
SA QN around the folding bands. As an example, 
a 3rd order GCF filter provides a EA QN rejection 
8dB higher than that guaranteed by a classical 3rd 
order comb filter. Throughout the paper we will use 
this optimal choice of the zero rotations where no 
otherwise specified. We invite the interested readers to 
refer to [9] for further details about the characteristics 
along with the performance of GCF filters. 

III. Partial Polyphase Decomposition of 
GCF Filters 

This section presents a non-recursive, partial 
polyphase implementation of GCF filters, suitable for 
decimation factors D that can be expressed as the p-th 
power-of-two, i.e. D = 2 P with p a suitable integer 
greater than zero. For conciseness, we only address 
the implementation of a 3rd-order GCF filter, but the 
considerations that follow, can be easily extended to 
higher order GCF filters with D as specified above. 

Let us focus on the optimal zero rotations shown in 
TableU and consider the z-transfer function Hgcf 3 ( z ) 

'For conciseness, we only deal with the positive semi-plane in 
the z-domain; however, the zero placement is specular for what 
concerns the zeros located in the lower semi-plane. 
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derived in (0. Due to the symmetry of the q p s in 
Table U the zeros belonging to Hgcf 3 { z ) can be 
collected in the following three z-transfer functions: 
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whereby (3% = z-e JQ , fa = z-e~ JC \ and a — \qi\2irf c . 
Notice that the first z-transfer function accounts for the 
zeros falling in the digital frequencies yj, k G Kk, 
which corresponds to the zeros of a classical 1st order 
comb filter. The other two z-transfer functions consider 
the rotated zeros. 

The z-transfer function Hgcf 3 {z) can be easily 
obtained by multiplying the previous three z-transfer 
function^]: 

(5) 

The impulse response of the 3rd order GCF filter in (|5]l 
has been derived in the Appendix. 

A non-recursive implementation of filter Hgcf 3 (z) 
can be obtained by expressing each rational function 
in © in a non recursive form. By doing so, the first 
polynomial ratio can be rewritten as follows: 

D-l log 2 (D)-l 

= £«-*= n O+^l «> 



1 



1-2!" 



i=0 i=0 

whereby last equality holds for D = 2 P . Upon using a 
similar reasoning, it is straightforward to observe that 
the following equality chain, which derives from ^ 
by imposing (3 — z ■ e T ^ a , holds as well: 



T lo g2(-D 
=0 



n 



1 



±j2 i a 



(7) 



2 Notice that, for conciseness, in the mathematical formulation that 
follows, we omit the constant term assuring unity gain at base-band. 
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Fig. 3. Architecture of the partial polyphase implementation of the decimation filter Hqq p 3 (z). The overall decimation factor D = Di -D2 
is split between the polyphase section decimating by D± = 2 p p +1 , and the non recursive section, decimating by D2 = 2 p_p p~ 1 , composed 
of p — p p — 1 decimation stages each one decimating by 2. Integer p p can take on any value in the set [—1, p — 1], 



By noting that 

= 1 + 2cos{2 l a)z- 2 ' + z- r+1 
after some algebra, (0 can be rewritten as follows: 



T log 2 (£>)-l 



l + z- 



Hgcf 3 (z)=UZ! 
■ (l + 2cos(2 t a)z- r +z- 2l+1 ) 

= YtS (D) - 1 [l + 
whereby 



whereby r = 1 + cos(aDi), and a = gi 1 7r -i . 

The polyphase decomposition of the z-transfer func- 
tion Hp(z) is defined as follows: 

H P (z) = Eklo 1 z- k E k (z Dl ) (13) 

whereby the functions Ek(z) are the polyphase com- 
ponents. Time-domain coefficients related to Ej.(z) are 
defined as follows: 



e k (n) = h P (D in + k), Vfc G [0, D x - 1] 



(14) 



(8) 



n = l + 2 cos (2*a) = 1 + 2 cos (g 2 ^ 
= 1 + 2 cos (q2 i+1 nf c ) , Vi = 0, . . .,log 2 (D) - 1 

Assume that the decimation factor Z? can be decom- 
posed as follows D — D± ■ D 2 , whereby £>i = 2 Pp+1 
and £>2 = 2 p ~ Pp ~ 1 . By doing so, © can be rewritten 
as follows: 



HgcfM = H P (z)-H N (z) (9) 



whereby 



Hp(z) = UZo 

Ip-] 
*=pp+i 



1- 



-3-2* 



and can be easily obtained by employing the recursive 
equation in (fTTT i. 

Some observations are in order. By choosing p p — 
p — 1, GCF filter Hgcf 3 (z) is fully realized in 
polyphase form, whereas for p p = — 1 the overall 
decimator is realized as the cascade of p non recursive 
decimation stages each one decimating by 2. Any other 
value of p p G [0,p — 2} yields a partial polyphase 
decomposition. 

The natural question that arises at this point con- 
cerns the practical implementation of the GCF filter 
Hgcf 3 (z)- In the following we derive an architecture 
for implementing GCF filters, while in the next section 
we present a mathematical framework for highlighting 
the sensitivity of the proposed architecture to the 



■ z 
(10) 

Remembering that p p — log 2 (_Di) — 1, it is straightfor- 
ward to observe that Hp(z) is the z-transfer function 
of a 3rd order GCF filter decimating by D\. The 
impulse response hp(n),Wn G [0, 3£>i — 3], of filter 
Hp(z) has been derived in Appendix: 



hp(n) 



— e +jo-n 



,-2 4 1 ,-2-2* \ 1 ,-3-2* approximation of its multipliers. The latter is needed 
for deducing useful hints at the basis of multiplier-less 
implementations of the proposed filters. 

By applying the commutative property employed in 
[5], it is possible to obtain the cascaded architecture 
shown in Fig. [3] The first polyphase decimation stage 
allows the reduction of the sampling rate by D\, 
thus reducing the operating rate of the subsequent 
n k 3 k 2 decimation stages belonging to Hm{z). Any stage of 

e -2j«fe3 ^ e jafc 2 ^ Xt ( fcl ) (H)h n (z) in Fig. [3]is constituted by a simple FIR filter 
fe 3 =o fc 2 =o fci=o operating at a different rate. Such an example, the i- 

th stage, with i G [0,p — p v — 2], is characterized 



The sequence Xt(n) is defined as follows: 

x t (n) = 8(n)-r5(n~D 1 )+r5(n~2D 1 )~8(n-'iD 1 ) 

(12) 



by the transfer function [l + r, + z~ 2 ) + z~ 3 ] 
operating at rate f s / • 2 l J (f s is the EA sampling 
frequency). 
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The frequency response related to Hn(z) can be Equation d22l can be rewritten as follows: 



evaluated by substituting z 



in ([Toll: 



■ [cos (3 • 2 u - 1 uj) + r u cos (2 u ~ 1 tjj)] 



(15) 



whereby uo = lirfd, and r u is defined as: 
r u = l + 2cos(2"a), Vu = p p + 1, . . . ,p - 1 (16) 

IV. Sensitivity Analysis 

This section deals with the analysis of the sensi- 
tivity of filter Hgcf 3 (z) to the approximation of its 
multipliers. In brief, the goal is to design the non- 
recursive architecture in Fig. [3] without multipliers, 
while guaranteeing the gain (8 dB based on the results 
shown in Table UJ) in terms of EA QN rejection with 
respect to a classical comb filter. 

Let us evaluate the sensitivity of the frequency 
response in (O with respect to its coefficients. Notice 
that there are two sets of coefficients: L — 3Z?i — 2 
multipliers (i.e. c n fc = hp(Din + k)) belong to the 
polyphase section Hp(z), andp— p p — 1 multipliers r u 
(shown in (TToT i) belong to the decimation filter Hn(z). 

First of all, notice that when the generic coefficient 
c n ,k (r u ) is approximated, its actual value can be 
expressed as c n , k = c n ,k + Ac„ tk (r u = r u + Ar u ), 
whereby Ac Uy k (Ar u ) is the approximation error. On 
the other hand, the approximations of coefficients c n> k 
and r u imply that Hgcf 3 be written as 

H G CF 3 (e ju ) = H GC F 3 (e ju ) + AH GC F 3 {e jul ) (17) 

The dependence of the frequency response 
Hgcf 3 (s juj ) on the approximation of its multipliers 
can be evaluated by differentiating ©: 

AH GC F 3 (en = H N (e^)AH P (e^)+H P (e^)AH N (e^ 

(18) 

whereby 

p-i 

,111-.. i , •-■ . 

-Ar„ 



dH N {e^) _ n r p ju>\ cos ( 2 " S 

dr u Ll N\^ ) cos (3.2"-iw)+r„-cos(2"- 1 w) 

(23) 

Upon substituting ( 1231 in (fT9l l. it is possible to 
obtain: 

H P {e 3U )AH N {e 3UJ ) = H P (e juJ )H N (e^)- 

^— \p— 1 cos(2 u ~ 1 ^j) Ar u 

' zJi4=p p + l cos(3-2"- 1 w)+r„-cos(2"- 1 w) = 

= H (p3 u \ VP^ 1 cos ( 2 "~ 1 ") Ar " 

J^GCF 3 \^ ) l^ u=Pp + i cos(3-2"- 1 w)+r u -cos(2"- 1 w) 

(24) 

Let us consider AH P (e 3UJ ). Given H P {e 3ul ) in (ED, 
it is straightforward to obtain the following relation: 

dH P (e juJ ) 



dc ntk 



By substituting the previous equation in d20t . it is 
possible to obtain: 

Dx-l VL/D 1 \ 

AH P (e 3UJ ) = J2 E Ac n , k e- MDin+k} 

fc=0 n=Q 

Upon multiplying and dividing by Hp(e 3UJ ), the func- 
tion H N (e ju )AH P (e juJ ) can be rewritten as: 



TT (f ju\ ^l- 1 ^/^iJ 



k=0 n=0 



(25) 



The actual frequency response HGCF 3 (e 3UJ ) in (fTTI i 
can be expressed as follows: 

Hgcf 3 {^) = H GC F 3 {e^) + AH GC F 3 {e^) = 



1 + Y P ~ 1 

J ^ Z-^i=Pp + l 



s(2 i_1 cj)Ar i 



=Pp + l cos(3-2* - ! u]+r, -008(2' 

(26) 

The effects of the approximation of the multipli- 
ers c„ fe and r u on the actual frequency response 
Hgcf 3 (c 1uj ) can be understood by analyzing the fre- 
quency behavior of the following error function: 



u=p p + l 



dr u 



(19) 
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H N (e^)AH P (e^) = H N {e 3 
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■e 
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ju(D in +k)] _l STP- 1 cos(2 i - 1 o;)An 
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dH P (e ]u 
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with 



Di-l [L/Di] 



—ju{Din+k) 



(20) 



(21) 



fc=0 n=0 



Let us evaluate the derivative of Hjs[(eJ u ) with respect 

to r u , \/u = p p + l,...,p-l: 

g-is-a™-^ ^ cos (-3 . 2™-^) + r m ■ cos (2™"^)] 

(22) 



'i=p P + l cos(3-2 i ~ 1 w)+r i -cos(2 i - 1 w) 
=P P H 

(27) 

which, to some extent, quantifies the distortion be- 
tween the desired HGCF 3 (e JUJ ) and the actual fre- 
quency response HGCF 3 (e 3u> )- 

Fig. [4] depicts the frequency behaviours of the error 
functions AHi (e juJ ) and AH 2 a (e JU ) noted in the 
last row of ( 1271) . for the following sample set of 
parameters: D = 32, D\ = 8, D 2 — 4, p p = 2, 
p = 5,v = 4,q=\qi\ = 0.79, and Ah P (n) = Ar^ = 
10" 4 , Vn,i. 

Fig. shows the behaviours of the frequency re- 
sponses Hp(e Jul ) and H^r{e 3U ) for the sample set of 
parameters noted in the respective label. 
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Fig. 4. Frequency behaviours in dB of the functions | AHi \ 
| Aff 2 ,3 (e J ") |. and | AH 2 ,4 (e jw ) | in (27) for D = 32, £>i 
8, D 2 = 4, p p = 2, p = 5, v = 4, g = 0.79, and Ah P (n) 
An = 10" 4 , Vn,i. 



Fig. 5. Frequency behaviours in dB of the functions \Hp 

andl_ff]v (eJ") I, respectively, in (13} and (BJ for D = 32, D\ ■. 

8, D 2 =i, Pp = 2, p = 5, u = A, q = 0.79. 



Some key observations are in order. Fig. |4] shows 
that the pass-band behaviour of the filter Hgcf 3 {& u ) 
is not affected by the approximation of its coefficients. 
Sensitivity of HacF 3 (e-'' J ) is very low for fd E [0, f c ], 
whereby f c — j-. This in turn suggests that the 
filter pass-band droop does not degrade by virtue of 
multipliers' approximations, and it is as low as the 
one guaranteed by filter Hgcf 3 (^ u )- Notice also that 
the sensitivity is very low around the digital frequency 
jj — fc, which defines the selectivity of the decimation 
filter [1]. 

Fig.sH and [5] show that the approximation of mul- 
tipliers hp(n) and r u affects the sensitivity of the 
frequency response Hgcf 3 {& iu ) m disjoint folding 
bands. Indeed, filters H P (e> u ) and H N (e^) place 
the respective zeros in different digital frequencies, as 
clearly highlighted in Fig. 

Fig. [4] also shows that the frequency response 
Hgcf 3 (&' u ') is very sensitive to coefficients' 
approximations specially in the folding bands 
[^-Zcs+Zc], Vfc G K k . This in turn 
suggests that particular care must be devoted to the 
approximation of multipliers embedded in both Hp(z) 
and Hn(z) in order to preserve the QN rejection 
performance around the folding bands. However, the 
same figure suggests that the approximation of the 
multipliers r„s belonging to H^(z), can be done 
independently from the approximations of coefficients 
hp(n) belonging to Hp(z). 

Sensitivity analysis derived above, allows to draw 
a general picture of the effects of the approximations 
of the coefficients on the actual frequency response 
Hgcf 3 (z)- Nevertheless, we deduced that the sensitiv- 
ity is very low in the pass-band [0, f c ] and around the 
frequency — f c . This in turn suggests that both pass- 



band droop and selectivity of GCF filters are preserved 
by the approximations of the multipliers. 

An important question is still open. We still need 
to quantify the extent of the effects of the coefficient 
approximations on the frequency response Hgcf 3 {z). 
In brief, the basic question we want to answer is as 
follows. What are the approximation errors Ar u , Vu = 
p p + l,...,p-l, and Ahp(n), Vn € [0,3A - 3], 
that we can tolerate on Hgcf 3 (z)1 The answer to 
this question is the focus of the next section. It is 
anticipated that whatever the condition on the maxi- 
mum approximation error tolerated on the frequency 
response Hgcf 3 {^ u ), it should be related to the 
behaviour of such a function around the folding bands, 
since both pass-band droop and selectivity of these 
filters are mainly unaffected by the approximation of 
the multipliers. 



V. Estimation of zero displacements due to 

COEFFICIENT APPROXIMATIONS 

Besides improving the selectivity on the frequency 
■fi—fa GCF filters provide improved SA QN rejection 
around the folding bands with respect to classical, 
equal-order comb filters [9]. However, coefficient ap- 
proximations can have detrimental effects on the zero 
locations in the z-plane, and, accordingly, can worsen 
SA QN rejection performance around the folding 
bands. As a consequence, it is useful to estimate the 
effects of coefficient approximations on the actual 
SA QN rejection performance guaranteed by filter 
HacF^eP^) by estimating the induced zero displace- 
ments. Next three sections address this topic by first 
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examining filters Hp(z) and Hpj(z) separately and 
then deducing some hints from the derived theoretical 
analysis. 

A. Displacements of zeros belonging to Hp(z) 

In this section, the focus is on the evaluation of the 
errors Az k on the locations of the zeros of Hp{z) 
due to the approximations of the coefficients hp(n) in 
the polyphase filter Hp(e^). First of all, notice that 
Hp{z) places its 3 • D% — 3 zeros in the following z 
locations: 



Zk = < 



i 

+j2*-(l5t+9/c) _ 
i27r(^-+g/c)_ 
±i*r(lfc-«/o) 



Vfc 


= 1,..., 






Vfc 


= 1,..., 


L^J - 


1 


Vfc 


= 1,..., 


L^J- 


1 


Vfc 


= 1,..., 


L^J- 


1 


Vfc 


= !,•••, 








in O 


can be 


ex 



pressed in a form emphasizing its zeros: 

L-l 

Hp{z)= [] (1-ZkZ- 1 ) 



fc=l 



whereby L — 3 • D\ — 2. Due to the approximation 
of the coefficients hp(n), the set of zeros becomes 
{z k = z k + Az fc , Vfc = 1, . . . , L - 1}: 



L-l 



ffp(^)=n( i - 



fc=l 



Zero displacements Az^ can be related to the coeffi- 
cient approximations Ahp as follows: 



L 



9zi 



A/i P (ry) 



(29) 



since there are L coefficients hp(n) with n € [0, L — 
1]. 

Upon noting that: 

dH P {z) dHp(z) 

U=. 



9zi 



dhp(rj) , _^ <9z 
it is straightforward to obtain: 

by employing ( [T3l l. and 

9i? P (z) 



dh P (ri) 



v 



L-l 



fe=i 



L-l 

n (i 

i=i,i=£k 



U=*< 



3 This is possible by virtue of the sensitivity analysis derived 
above: zeros of both Hp(z) and Hj\r(z) affects different folding 
bands. In other words, each pair of conjugate complex zeros affects 
the behaviour of the frequency response in only one folding band. 




Fig. 6. Maximum zero displacement over all the set of coeffi- 
cients hp(n) as a function of the approximation error Ahp. The 
results assume the same error Ahp for all the coefficients hp(n). 
Parameters are as follows: D2 = 4, and v = 4. 



by deriving ( 1281 1 with respect to z. Finally, after some 
algebra d29l ) can be rewritten as follows: 



(28) A Zi 



L 



zpAhpin) 



^{E^nf^W 1 -*'* -1 )}, 

(30) 

which is the displacement of the i-th zero Z{ due to the 
approximations of all coefficients hp(n) belonging to 
the polyphase filter decimating by D\. 

Fig. [6] shows the maximum Azi over all the set of 
zeros indexed by i = 1 , . . . , L, as a function of the 
approximation error Ahp(rf) = Ahp, assumed to be 
the same for all coefficients, for various values of D\, 
D-2 = 4, and v = 4. Curves shown in the figure can 
be considered as the worst case zero displacement due 
to the approximations of coefficients hp(n). 

B. Displacements of zeros belonging to Hn(z) 

In this section, the focus is on the evaluation of 
the displacements Azk of the zeros in Hjy(z) due to 
the approximation of the coefficients r u belonging to 
Hm{z). Error Azi on the i-th zero can be related to 
the approximation errors of multipliers r u s as follows: 



A Zl 

Upon noting that: 

dH N (z) 



p-i 



dzi 



^ dr, 

u=p p + l 



Ar„ 



(31) 



dr u 



dH N {z) 
dz 



dz l 
dr u 



after some algebra on ( flOt . it is possible to obtain: 



dH N (z) 



-2-T 



1 lk=p p -\-\ J kjby l I 



+r k (z- 2 " + z- 2 - 2 "^ + z- 3 ' 2 " I 



(32) 
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Fig. 7. Maximum zero displacement over all the set of coefficients 
r u as a function of the approximation error Ar. Curves assume the 
same error Ar for all the multipliers r u . Parameters are as follows: 
D2 = 4, and v = 4. 



By deriving ( TTOb with respect to z, it is possible to 
write: 



dH N (z) 

-2-2 k z~ 2 - 2 ' 



-3-2 k z 



-3-2" 



p-1 



i + n(z 



')+z-^}} 



(33 = f 

Finally, ( f3TT > can be evaluated by substituting the ratio 
between d32l > and (l33l in place of 

Fig. [7] shows the maximum Azi over all the set 
of zeros belonging to Hn(z), as a function of the 
approximation error Ar u — Ar, assumed to be the 
same for all multipliers, for various values of D\, 
D 2 = 4 and v = 4. 

A quick comparison between the results shown in 
Fig.s [6] and [7] reveals that the maximum zero dis- 
placement of the zeros belonging to Hm{z) is much 
smaller than the one experienced by zeros belonging 
to Hp(n). There are at least two basic reasons for 
such a behaviour. First, the number of multipliers 
r u s is very small with respect to the number of 
coefficients hp (n) of the polyphase section. Secondly, 
any error Ar u slightly rotates the zeros belonging 
to the u-th decimation cell in H^(z) leaving them 
on the unit circle. This follows from the z-transfer 
function [l + r u + z~ 2 ) + z~ 3 ] of the u-th dec- 
imation stage. On the other hand, approximation errors 
Ahp(n) can also move the zeros of Hp(z) outside 
the unit circle in the z-plane worsening the EA QN 
rejection performance of the decimation filter. 

The previous analysis suggests that the most crit- 
ical filter in the partial polyphase architecture is the 
polyphase section Hp(eP u ). In order to deduce the 
maximum tolerable approximation error over the set of 
coefficients Ahp(n), the behaviour of filter Hp{e> w ) 




Fig. 8. Behaviour of the difference AP qn = P qn \ dg — Pqn\ dB 
as a function of the decimation factor D\ and for various values of 
the approximation errors Ahp shown in the legend. 



around the folding bands should be further investi- 
gated. This is the topic addressed in the next section. 

C. Design Considerations 

The EA QN power falling inside the folding bands 

[jJT ~fc, 7JT + /J ' Vfc = ' ' ■ ' LlH can be define <fl 
as [1]: 



fe=l 



-+fc 



\Hp(en\ 2 S B (fd)df d (34) 



where Spifd), the power spectral density of the 
SA QN, can be expressed as Spifd) = S e (fd) • 
[2 sin(TT fd)] 2B ■ In the previous relation S e (fd) — j^j- 
is the power spectral density of the sampled noise 
under the hypothesis of representing the QN as a white 
noise [1], A is the quantization level of the quantizer 
contained in the EA modulator [1], and f s is the EA 
sampling rate. 

In order to quantify the EA QN rejection perfor- 
mance of filter Hp(z) embedding the approximated 
multipliers hp(n) — hp(n) + Ahp(n), with respect 
to a classical 3rd order comb filter, the following 
performance metric, AP qn , can be evaluated: 



AP, 



qn 



D Jc 



Ud)\ 2 S B Ud)df d 



(35) 



The behaviour of AP gll -[dB] as a function of D\ 
is shown in Fig. [8] for various values of Ahp, i.e., 
the approximation error which is assumed to be the 
same for all the coefficients hp(n). Results in Fig. [8] 
show that EA QN rejection improvements can still be 

4 Note that j34t is valid only if the Noise Transfer Function (NTF) 
of the modulator is maximally flat, i.e., it does not contain stabilizing 
poles. In higher order modulators this requires multi-bit feedback 
structures. 
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Fig. 9. Modulo of the frequency response of filter H^cp^^^) 
employing approximated coefficients hp(n) = hp(n) + Ahp(n) 
(continuous curve), and frequency response of filter i?GCF 3 (e J ") 
embedding the real coefficients hp(n), for D = D± = 64, u = 4, 
and Ahp(n) = Ahp = 10 — 3 , Vra. Lower subplot shows the 
behaviours of such functions around five folding bands. 

achieved upon approximating each coefficient hp(n) 
within an error less or equal to 10~ 3 . 

This in turn suggests that it is possible to optimize 
each multiplier hp(ri) in the polyphase filter Hp(z) 
by approximating it as a power-of-2 (P02) coefficient 
with an approximation error of 10~ 3 without affecting 
the EA QN rejection performance of the polyphase 
filter around the folding bands for any D\ > 32. 

Fig. [9] shows the behaviours of both the frequency 
response of filter HQ CF3 {e^ u ) employing approxi- 
mated coefficients hp(n) — hp(n) + A.hp(n), and the 
frequency response of filter Hggf 3 {^ u ) embedding 
the real coefficients hp(n), for D = D-y = 64, 
v = 4, and Ah P (n) = Ah P = 1CT 3 , Vn. This is the 
most critical case in which a full polyphase architec- 
ture is employed. The lower subplot shows the local 
behaviour of both frequency responses around some 
folding bands. Notice that both curves are superim- 
posed, even though H GCF;s (e : * u> ) employs coefficients 
approximated with an error equal to 10~ 3 . 

The approximation of the multipliers hp(n) with 
P02 coefficients is beyond the scope of this work. 
Nevertheless, many excellent works have been pro- 
posed in literature for obtaining the best P02 coef- 
ficient approximation within a predefined error. We 
invite the interested readers to refer to papers [17]- 
[19]. 

Let us summarize the main design considerations 
deduced from the proposed sensitivity analysis of 
partial polyphase GCF filters. 

• Both pass-band droop and selectivity performance 
of GCF filters are mainly unaffected by the 
approximation of the multipliers embedded in 
the decimation filter. Practically speaking, this 



means that such performance are the one already 
deduced in the companion paper [9]. 

• The most critical section in the partial polyphase 
decomposition is the polyphase filter Hp(z). 

• Upon approximating each multiplier of the 
polyphase section with an approximation error 
less or equal to 10 -3 , it is possible to preserve the 
SA QN rejection performance of GCF filters with 
respect to classical, equal-order comb filters, for 
a wide range of decimation factors Di as noted 
in the abscissa of Fig. [8] 

VI. Conclusions 

This paper focused on the design of generalized 
comb filters by proposing a novel partial polyphase 
architecture with the aim to reduce the data rate after 
the EA A/D conversion. We proposed a mathematical 
framework in order to analyze both the sensitivity of 
the frequency response and the displacements of the 
zeros in the filter transfer functions due to the quan- 
tization of the multipliers embedded in the proposed 
filters. 

We also derived the impulse response of a 3rd order 
sample GCF filter, that we used as a reference scheme 
throughout the paper. 
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Appendix 

In this Appendix, we derive the impulse response 
hacF 3 (n), Vn = 0, . . . , 3D — 3, of a 3rd-order GCF 
filter decimating by D. The proof relies on repeated 
applications of the inverse z-transform on the product 
of two analytical functions X\(z) and X 2 (z) related 
to two discrete-time sequences, x\(n) and x 2 (n): 



Z- 1 [X 1 {z)X 2 (z)}=x 1 (n)-kx 2 {n) 
= Tl=-oo xi{k)x 2 {n - k). 



(36) 

First of all, consider the transfer function in ©, and 
define as X t (z) the numerator z -polynomial: 

X t (z) = (l - z- D e^ aD ) (1 - z- D e-i aD ) (l - z' D ) 



With the setup above, Hgcf 3 {z) m © can be 
rewritten as follows: 

HgcfM = X t {z) ■ Yx{z) ■ Y 2 (z) ■ Y 3 (z) 

Upon applying d36l l to the z-function Wi(z) = 
X t (z)Yi(z), it is possible to obtain: 

wi(n) = ^2 x t {kt)yi{n - ki) (38) 

k\— — oo 

Applying ( f36l > to the z-function W 2 (z) — Wi(z)Y 2 (z), 
and employing Q81l . it is possible to obtain: 

w z( n ) =Efc 2 =-oo W l( fc 2)y2(™-fc2) = 

= Efe^-txjEfe^-oo^OyiC^ - h)y 2 (n - k 2 ) 

(39) 

Finally, applying (l36l l to the z-function Hgcf 3 {z) = 
W 3 (z) = W 2 (z)Y 3 (z), and employing (l39l . it is 
possible to obtain: 

hacF 3 {n) = E w ^h)v^{ n - k 3) = 

&3 — — OO 

+oo +oo +oo 

= E E E M*i)yi(fc-*i)- 

fca= — oo k2 — — oo fc i — — oc) 

■V2{h - k 2 )y 3 (n - k 3 )} (40) 

Upon substituting the respective expressions of the 
sequences yi(n), Vi — 1, . . . , 3 in d38l ). it is possible 
to write: 

+oo +oc +oo 

hGCF 3 ( n ) = E E E [ x *( fc l) u ( fc 2 - fcl)- 

A,'3 — — oo k2 — — oo k\— — oo 

. e -M*3-k2) u ( k3 _ k 2 )e +3a{n - k3) u{n - k 3 ) 

By exploiting the definitions of the unitary-step se- 
quences, it is possible to observe the following rela- 
tions: 

1 Vfc 2 > h 



[l — rz 



-D i rz -2D _ Z -3D1 



whereby r = 1 + 2 cos(aD). The discrete-time, causal 
sequence with z -transfer function X t (z) can be written 
as follows: 

x t (n) = 8(n) -rS(n-D) + r8(n - 2D) - 5{n - 3D) 

(37) 

whereby S(n) is the discrete-time unit impulse cen- 
tered in n = 0. 

Let us define the following pairs of transfer func- 
tions along with the respective discrete-time sequences 
[16]: 

l 



u(k 2 — ki) 
u(k 3 - k 2 ) 
u(n - k 3 ) 



Vfc 2 < ki 

1 Vfc 3 > k 2 
Vfc 3 < k 2 



1 Vn > k 3 
Vn < k 3 

Based on these observations, we can reduce the 
upper limits of the summations in fill as follows: 

n k 3 k 2 

h G CF 3 (n)= E E E M fc i)' 



fcg — — OO k2 — — OO fci— — OO 

. e -3<x(k3-&2) e +jot(n-k 3 ) 



(42) 



Yi(z) 
Y 2 (z) 
Y 3 (z) 



2/i (n) = u[n) 
y2(n) — e~ ja u(n) 
y 3 (n) = e+3 a u{n) 



whereby u(n) is the discrete-time unitary-step se- 
quence. 



By observing that the sequence Xt(ki) is causal 
(see (|37|>). i.e., x t (ki) = 0, Vfci < 0, it is possible 
to obtain: 

n k 3 k 2 

h G CF 3 {n) = e +jan e ~ 2lak3 E ^ E Xt ^ 

k 3 =0 fe 2 =0 fei=0 
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Notice that the impulse response hGCF 3 (n) general- 
izes the one of the classical 3rd order comb decimation 
filter, and, it is composed by 3D — 2 coefficients 
over the time interval ranging from to 3D — 3. 
Equ. j43b can also be used as an on-line algorithm 
for generating the coefficients of the GCF impulse 
response hccF 3 (n) by simply solving the three nested 
summations for each n = 0, . . . , 3D — 3. 

As a note aside, notice that by imposing a = 
in d43l l. it is possible to obtain the impulse response 
of a classical 3rd-order comb filter. 
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